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In this work, a convection-driven Rijke-Zhao thermo-acoustic-piezo system is designed and experimen¬ 
tally tested to demonstrate its potential for harvesting thermal energy. For this, a nonlinear theoretical 
model is developed to simulate the energy conversion process, i.e. heat-to-sound-to-electricity in the 
present system. Unlike the conventional conduction-driven thermoacoustic converters, our present 
system involves no heat exchangers and stacks. As a heat source is placed in a Rijke-Zhao tube with 
two bifurcating daughter branches, self-sustained thermoacoustic oscillations are generated. The resulting 
acoustic fields in the bifurcating branches are dramatically different. One branch is associated with ‘hot’ 
oscillations. However the other is with ‘cold’ oscillations at ambient temperature, which enable a piezo¬ 
electric generator being implemented to the end of the branch. In order to measure the acoustic fields 
in the bifurcating branches, two arrays of thermocouples and microphones are used. The maximum sound 
pressure level is around 139 dB. The output electric power and acoustical energy conversion efficiency are 
measured and compared with that from a similar but a conduction-driven thermo-acoustic-piezo system. 
It is found that 60% more power is generated. And the energy conversion efficiency is increased by 105%. 
These experimental results confirm that the developed Rijke-Zhao thermo-acoustic-piezo system is an 
invaluable tool in designing a simple, low-cost, energy-efficient thermoacoustic system. 

© 2012 Elsevier Ltd. All rights reserved. 


1. Introduction 

Greenhouse gas emission and fossil fuel exhaustion are increas¬ 
ingly urgent issues. However, industry waste heat generated by 
burning fossil fuel is released to the atmosphere. This is not only 
harmful to the environment, but also energy wasting. To recover 
waste heat, there is a need for energy conversion techniques to 
convert heat into electricity or mechanical work, especially for heat 
at low temperature. Thermo-acoustic energy conversion technol¬ 
ogy [1,2] may be used for energy harvesting. Generally, thermoa¬ 
coustic energy conversion systems are associated with two 
distinct types of oscillations: (i) convection-driven, or (ii) conduc¬ 
tion-driven. Rijke tube [1-6] is a typical example of convection-dri¬ 
ven thermoacoustic systems. It is a simple open ended vertical 
tube with a confined heat source in its lower half. And a net flow 
of gas through the tube is generated due to the natural convection, 
which is essential in producing self-sustained oscillations. Thus the 
Rijke oscillations are convection-driven. Sondhauss tube [4,5] is a 
classical device of producing conduction-driven oscillations. It is 
a pipe having one closed and one open end. When heat is added 
and transferred through the wall of the closed end via conduction, 
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Sondhauss oscillations occur. There is no net flow of gas through 
the tube and so the Sondhauss oscillations are conduction-driven 
in the presence of mean temperature gradients [7[. 

To elucidate thermoacoustic energy conversion process [8,9], as 
occurred in Rijke- and Sondhauss-tube as shown in Fig. 1, consid¬ 
erable research has been carried out. The first insight into thermoa¬ 
coustic oscillations was provided by Rayleigh [10] over a century 
ago. He stated that intensive thermoacoustic oscillations would 
occur when heat was added in phase with pressure, but would 
not occur if heat was added out of phase with pressure. Culick [11] 
derived a general form of the Rayleigh’s criterion, which is applica¬ 
ble to both linear and nonlinear thermoacoustic oscillations in any 
shaped chamber. He showed that an oscillating heat source 
functioned like an oscillating piston. The same conclusion was 
drawn by Chu and Ying [12] and Rott [8], Recently, Bisio and 
Rubatto [6] reviewed the heat-generated acoustic oscillations 
in Sondhauss and Rijke tubes and discussed their possible 
applications. 

Besides offering insight into the energy exchange that gives rise 
to and limits the growth of intensive oscillations, Rayleigh criterion 
also suggests ways in which thermoacoustic oscillations can be 
maximized or minimimized, depending on applications in practice. 
One application in aero-engines and gas turbines is to minimize con¬ 
vection-driven thermoacoustic oscillations to stabilize combustion 
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Nomenclature 



a 

angle between the bifurcating branches (°) 

Pm 

surface density of the piezo diaphragm (kg/m 3 ) 

A f 

flame surface area (m 2 ) 


damping of the piezo diaphragm 

Ai, A 2 , A 3i A 4 cross-sectional areas of different regions along the 

t 

time (s) 


tube (m 2 ) 

T 

temperature (°C) 

Amn 

coefficient of the green function 

m 

tension of the piezo diaphragm 

C 

speed of sound in air (m/s) 

Cgt Tft, Tj, 

Tf, t time delay between incident and reflected 

Cm 

sound speed in the membrane (m/s) 


waves (s) 

c p 

specific heat capacity at constant pressure 

0 

estimated phase between p'(zpt) and Q(t) (°) 

El 

electrical power output (mW) 

U 

mean flow speed (m/s) 

>h 

acoustical energy conversion efficiency 

Ug 

oncoming flow velocity at flame position (m/s) 

mm 

incident and reflected waves in region 1 (Pa) 

U, V 

coefficient matrices 

g 

green function 

V 

volume flow rate (m 3 /s) 

G (z, r, t) 

flame front position 

Vrms 

root mean square of the output voltage (V) 

mj(t ) 

incident and reflected waves in region 2 (Pa) 

Vf 

laminar flame speed (m/s) 

H(co) 

transfer function matrix of the piezo transducer 

V(co) 

voltage generated by the piezo diaphragm (V) 

Jm 

Bessel function of order m of the first kind 

z 

axial position (m) 

Km 

stiffness of the piezo diaphragm (kN/m) 


, z d4 upstream and downstream open ends axial posi- 

K 

wave number 


tions (m) 

Kmn 

the nth root of Bessel functions of order m of the first 

Z u (co) , Z 3 ((u),Z 4 (m) reflection coefficients at the open ends 


kind 

z(oj) 

deflection of the piezo diaphragm (m) 

/(f). m(t) 

incident and reflected waves in region 3 (Pa) 



L(m) 

Fourier transform of incident wave L,(t) 



m m 

mass of the piezoelectric diaphragm (g) 

d 

downstream of the flame section 

M 2 

Mach number at region 2 

u 

upstream of the flame section 

n f 

the unit vector of the direction normal to the flame sur- 

i 

ith region along the tube 


face 

m 

piezo diaphragm 

CO 

radian frequency (rad/s) 

rms 

root mean square 

pM 

acoustic pressure (Pa) 



App(m) 

pressure difference across the piezo diaphragm (Pa) 


Q 

unsteady heat release rate (kj/s) 

- 


q(m) 

piezo-diaphragm frequency-dependent charge 

A 

Fourier transform 

m s(t) 

incident and reflected waves in region 3 

' 

fluctuation part 

Ri(t),R|(m) reflected waves in time- or frequency-domain 



R, 

load resistance (£2) 



PM 

air density fluctuation at ith region (kg/m 3 ) 




systems [13-16]. These have been well-reviewed by McManus et al. 
[17] and Dowling and Morgans [18]. The other application is to max¬ 
imize the thermo-acoustic oscillations to achieve energy conver¬ 
sion. Thermo-acoustic engines [19,20] are generally associated 
with conduction-driven oscillations. Such engines consist of an 
acoustic resonator, heat exchangers, a stack (or regenerator). And 
they can be divided into two varieties; one is standing-wave engine 
[21-23] and the other is travelling-wave one [24-26]. Both classes 
of engines have under vigorous development, after Ceperley [20] 
realized that Stirling engines are of the travelling-wave class. 

Thermoacoustic engines are more attractive than those tradi¬ 
tional Stirling engines involving pistons, due to the fact that they 
lack moving parts, require little maintenance and last a long time 
[27]. Different designs of thermoacoustic energy conversion system 
and various forms of input energy have been tested to generate 
large-amplitude or high-frequency sound waves [28,29]. Miniature 
thermoacoustic engines were also developed [30,31]. Optimization 
of thermoacoustic system performance was studied by Babaei and 
Siddiqui[32], Thermoacoustic converters powered by engine waste 
heat [33], or solar power [34,35] have been developed and tested. 
Recent developments [36] and application of thermoacoustic 
system for energy harvesting was experimentally investigated by 
Smoker et al. [23]. A conduction-driven thermo-acoustic-piezo sys¬ 
tem was designed. A porous stack and two heat exchangers were 
used to generate a steep temperature gradient to generate self-sus¬ 
tained thermoacoustic oscillations. The resulting oscillations were 
then used to excite a piezolelectric diaphragm placed at the end 


of the resonator. The output power is found to be very low (approx¬ 
imately 0.128 mW). 

In this work, a convection-driven Rijke-Zhao thermo-acoustic- 
piezo system is designed and experimentally tested to demon¬ 
strate the possible application for energy harvesting. For this, a 
nonlinear model is firstly developed to simulate the energy conver¬ 
sion process in the Rijke-Zhao thermo-aoustic-Piezo system. This 
is described in Section 2. The experimental setup and measure¬ 
ment of the Rijke-Zhao energy conversion system is then de¬ 
scribed in Section 3. Two arrays of microphones, thermocouples, 
an infrared camera and a constant temperature hot-wire anemom- 
etry are employed to monitor and analyze the acoustics fields. 
Finally, in Section 4, the electricity generated from the our present 
system is measured. Comparison is then made between our results 
and that of a similar but conduction-driven thermo-acoustic-piezo 
system as proposed by Smoker et al. [23], 


2. Model of Rijke-Zhao thermo-acoustic-piezo energy 
conversion system 

To investigate the feasibility of the harvesting thermal energy, a 
nonlinear thermo-acoustic model of Rijke-Zhao tube with a piezo¬ 
electric generator implemented is developed, as shown in Fig. 2. 
The model captures three physical processes; the generation of 
acoustic waves by a premixed laminar flame, the propagation of 
pressure waves [37,38], and the dynamic response of a piezoeletric 
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diaphragm to acoustic oscillations. It is worth noting that the pres¬ 
ent system is different from conventional Rijke tube due to the un¬ 
ique structure and geometric configuration and the resulting flow 
patterns in the bifurcating branches [39]. It will be shown later that 
the present system can produce both ‘hot’ and ‘cold’ thermoacou¬ 
stic oscillations, which enables the energy harvesting achievable 
by using piezoelectric generator. However, conventional Rijke tube 
produces only ‘hot’ oscillations. 


2.1. Model for pressure wave generation by unsteady heat release 


The heat source in our experiment is a premixed laminar flame. 
The presence of the heat source separate the Rijke-Zhao combus¬ 
tor into 4 regions, as shown in Fig. 2. The flame is assumed to be 
acoustically compact and so it is modelled as a thin sheet. The flow 
fluctuations upstream and downstream of the flame can be mod¬ 
elled as being due to plane acoustic waves travelling in opposite 
directions. Thus the equations for the flow fluctuations upstream 
and downstream of the flame are expressed in terms of acoustic 
waves R, and L, as given 


Pi(z,t)=p + Ri t- 


C i + Ui 


+ Ljl t + 


Ui{Z,t) =Ui + i 

, , l 

ft(Z, t) = Pi + =2 


R,[t- 
Rift — 


Cj + Ui 
2 


Ci + Ui 


Ci - Ui 
— Li ^ t + 

-u( t 


Ci - U 

z 


Ci - u, 


(la) 

(lb) 

(lc) 



where p denotes the instantaneous pressure, u is the flow velocity, 
p is the density, c the speed of sound and subscript i denotes the ith 
zone. Li denotes the incident wave g(t), j(t), m(t) or s(t), and Rj 
denotes the reflected one/[t), h(t), r(t) or I(t ) in the different regions. 
The instantaneous value consists of a mean value denoted by an 
overbar, and a perturbation part denoted by a prime. R,(t) = 
Ri(a>)e> mt and Lft) = Lj(a>)e> M denote the reflected and incident 
wave respectively. Here flow disturbances are assumed to be 
time-dependent e’ mt . 

The momentum and energy equations across the flame hold 
as 


1 

P 2 -Pi+p 2 u 2 (u 2 -u 1 )+-c D py i =0 ( 2 ) 

( p 2 u 2 -PiUi)+^p,Ui (ui -u])= Q- {t) (3) 

The term on right hand side of Eq. (3) describes the heat release 
from a given heat source. To capture the heat release dynamics, a 
nonlinear flame model is developed. The flame model describes a 
premixed laminar flame, burning in the wake of a metal wire gauze 
anchored to a Bunsen burner in Rijke-Zhao tube. The fluctuations of 
the oncoming flow velocity at the gauze result in flame area varia¬ 
tions and further lead to unsteady heat release, which is an efficient 
sound source. The instantaneous total rate of heat release Q(t) = 
Q + Q'(t) has been shown [14] to depends on the flame surface area 
Af(t) — Af + Aj as 

Q + Q'(t) A f + /\' f (t-x v ) 

Q. Aj 

where t v represents a time delay between Q(t) and Aft), which de¬ 
pends on the flame surface front. 

The conical flame surface is assumed to propagate itself in the 
direction of its normal %= VG/|VG|, into the unburnt fluid at a 
constant laminar burning velocity Vf relative to the unburnt gas. 
Here the flame surface axial position at radius r is given by 
G(z,r,t) = z- {(r,t) = 0. z is the axial coordinate and t is time. As 
shown in Fig. 2, the surface G = 0 moves in the direction of its 
normal with speed Uf— Vpf, where u f= (u, v) is the unburned gas 
velocity, and n f points upstream. Mathematically this statement 
is equivalent to dG/dt = 0: 

dG dG r 

dt = at + (U/ “ V/D/) ' VG (5) 


After substituting n/= VG/|VG| and G(z,r,t) = z - £(r,t), the familiar 
G-equation is obtained, which describe the location of the flame 
surface [14]. Substituting for G(x,r,t) = z - £(r,t) gives 


di 

— + v — 

dt dr 


+ V f 



= U 


( 6 ) 


Since the fuel-air ratio is so low for the laminar flame that the over¬ 
all expansion is negligible. Thus the density change across the sur¬ 
face G(z,r,t ) = 0 is assumed to be negligible. The particle velocity in 
Eq. (6) is given by that in the oncoming flow (u G ,0). Eq. (6) then 
simplifies:ce:display> 


d£ 

dt 


= u c -V f 



( 7 ) 


To solve Eq. (7) so as to determine the flame front location i[r,t) for 
a specified oncoming velocity, numerical method such as Runge- 
Kutta is needed. For simplicity the oncoming velocity is assumed 
to be independent of radius, i.e. u c (t) as given as 


u G = u c + 


1 

Pi Cl 


(f(f)-g(t)) 


( 8 ) 
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Once £(r,t) is known, the instantaneous flame surface area Aj(t) can 
be readily calculated, since 

w-M-FW* <9> 

Thus linearizing the momentum and energy Eqs. (2) and (3) gives 
the time evolution of the flow fluctuations f[t), g(t), h(t ) and j(t) 
generated by the unsteady heat release Q'(t). 

At the junction section of the bifurcating branches, mass and 
energy conservation hold as 

f, [ pdV = p 2 u 2 A 2 - p 3 u 3 A 3 + p 4 U 4 A 4 = 0 (10a) 

OL JdV 


Ft L [W +pc ^) Av = p^(l u §+ c ^) 

+ p 4 u 4 A 4 (^uj + c p T 4 ) - A 2 p 2 u 2 Q u\ + c p T 2 j =0 (10b) 


where A u A 2 , A 3 and A 4 denote the cross-sectional area of the com¬ 
bustor in different regions. Linearizing the mass and energy Eqs. 
(10a) and (10b) and expressing the flow fluctuations in terms of 
r(t), s(t), /(t) and m(t) gives the time evolution of the acoustic waves 
propagated in bifurcating upper branches. However, for conve¬ 
nience, pressure continuity can be assumed to hold at the junction 
section, i.e. 

— Zfi - jcu——— , 

h(co)e E 2 ( ,+M 2 )+j(cu)e £ 2 n-*» 2 > = r(co) + s(co) = l(co) + m(co) (11) 

where overhat denotes Fourier transform of the corresponding 
parameter in time domain. 

The combustor boundaries at the ends of the mother tube and 
the bifurcating branches are modelled using pressure reflection 
coefficients Z U ,Z 3 and Z 4 . In frequency domain, they correspond 
to Z u (co), Z 3 (a>) and Z 4 (co) respectively [14]. These reflection 
coefficients relate the reflected wave to the incident wave at the 
boundary, as 

gj.-2„ (0 ,)e--., 

= -Z 4 (co)e~ im ‘ ( 12 ) 

l(co) 


For convenience, the value -1.0 can be used at the open upstream 
and downstream end, i.e. Z U = Z 4 = -1.0. To make our analysis 
more general, the pressure reflection coefficient at the end with pie¬ 
zoelectric diaphragm attached is assumed to be Z 3 . 

Thus the end boundary conditions and the flow conservation 
equations across the flame and the junction section of the bifurcat¬ 
ing branches provide enough information to solve for each of the 
eight wave strengths in Fig. 2. The resulting matrix equation is 
given as 


(h(t)\ 


/ h(t - T h ) \ 


/ 0 \ 

g(t) 


g(t~ tg) 


0 

m 

= V 

j(t ~ *j) 

+ 

0 

r(t) 


r(t - T r ) 


0 

v m J 


V Kt-ti) ^ 




where U and V are coefficient matrices as given in Appendix. The 
time delays are given as 


Tu _ z d t — 2 Zu _ 

h c 2 (l+M 2 )’ * c,(l-A*f) 


r. -_ 

■ h — 


c 2 (1-M 2 ) 


= 2[Zfl— Zdj/CQS a and _ 2(Z J4 -Z d )/C0S« M the mean Mach num b er i n 
c 3 (1-MZ) 1 c 4 (l-Mj) 1 

ith region of the Rijke-Zhao combustor. 

Now the time evolution of the flow perturbations from specified 
initial conditions can be determined by numerically solving 


Eqs. (7), (9) and (13) using a fourth-order Runge-Kutta method. 
In our calculations, ^(r,f) was evaluated at 50 locations across the 
duct radius R = 0.025 m, the length of the mother tube is 0.3 m, 
the bifurcating branches are 0.5 and 0.9 m respectively and the 
flame is placed at 0.13 m away from the bottom open end. 

Fig. 3a shows that the time evolution of the heat-driven acous¬ 
tic oscillation. It can be seen that initial flow disturbance is quickly 
grown into limit cycle. The resulting thermoacoustic oscillation is 
at around 170 Hz with a sound pressure level 138 dB. Harmonics 
are also observed, as shown in Fig. 3b. Fig. 3c shows the insightful 
information about the nonlinear periodic oscillation in the phase- 
space diagram between pressure and velocity. Fig. 3d illustrates 
the variation of the phase difference 0{co 4 ,t) between the pressure 
p'(t) and heat release Q[(t) at frequency co with time, which is 
defined as 


8(co, t ) = arccos 



p'(co, 0)Q'(o>, 0)d0 

^fT /L ° 0)d0 O' 2 (co, 0)d0 y 

(14) 


It can be seen that when initial pressure perturbation grows into a 
limit cycle, unsteady heat release is in phase with the pressure fluc¬ 
tuation, i.e. \0(a> 4 , t)| r= 65° < 90°. This phase information is consis¬ 
tent with Rayleigh’s criterion [18], verifies our numerical model 
and reveals the generation mechanism of heat-driven acoustic 
oscillations. 

The time evolution of the heat-driven velocity oscillation is 
shown in Fig. 9. By comparing the estimation with the measure¬ 
ment, good agreement is observed. This validate our present 
heat-driven acoustic wave model. 


2.2. Model for electricity generation by deformed piezoelectric 
diaphragm 

Piezoelectric materials produce a voltage in response to an 
applied force because of charge redistribution caused by deforma¬ 
tion. Similarly, a change in dimensions can be induced by the appli¬ 
cation of a voltage to a piezoelectric material. In our application, 
the acoustic pressure oscillations is applied to a circular piezoelec¬ 
tric diaphragm, which is deformed and leads to voltage output. 
Following the work [23,40], the coupling between the acoustic 
and electrical fields from the a given point on the piezo-diaphragm 
is formulated into a two-port model in frequency domain as 

f z(co) \ = ('H u (co) H u (co)\ /Ap p (m)\ 

\Q(co)J \H 2 j(co) H 22 (co)J [ V(co) J (15) 

H(O) 

where z(a>).q(co),Ap p (a>) and V(co) denotes the piezo-diaphragm 
freqeuncy-dependent center point deflection, charge, pressure dif¬ 
ference and voltage. The transfer function matrix H(co) consists of 
the elements of the admittance of the piezo-diaphragm. Here 
Hu(a>) and H 22 (oj) are the acoustical and the electrical admit¬ 
tances, and H 2 i(cb) and H 12 (ct>) denote the coupled electro-acoustic 
admittances. Detailed information about how to estimate the trans¬ 
fer function H(co) can be found in the work of [23]. Manipulating 
Eq. (15) gives the transfer function between voltage output and 
pressure fluctuation as 

V(co) H 21 (co) — Hu (a>) 1 f z(to) q(m) 1 

a Pp(®) H 12 (co) -H 22 (ffl) H u ((d)-H 22 ((o) L a Pp(®) App(co)J 

The displacement of the thin circular piezoelectric diaphragm 
z(r, 0, t) = z{r , 6) exp(jmt) results from the pressure difference be¬ 
tween the two sides of the diaphragm, i.e. A p p (r.0,t) = 
App(r) exp(jmt), as shown, 




























D. Zhao/Energy Conversion and Management 66 (2013) 87-97 


91 


(a) 


a. 


5 


0 


-5 

0 0.5 1 1.5 2 


1 limit cycle 



i 

i 


Time (s) 


(b) 


150 


CQ 


100 


50 


o- 

G 

3 

O 


100 200 300 400 

Frequency (Hz) 





Fig. 3. Onset of Rijke-Zhao thermoacoustic oscillation: (a) Time evolution of pressure oscillation, (b) frequency spectra, (c) phase diagram between u G {t)/u c and p(Z/, t)/p and 
(d) phase diagram between Q(t)/Q and p(z/,t)/p. 


dt 2 + m m dt 



z(r, 8)e> mt 


£m 

Pm 


\d_ 
7 dr 




J _o7_ 

+ r 2 W 


(*M)) 


gjwt 


_ gjmt 

Pm 


(17) 


Eq. (17) is an inhomogeneous equation describing the dia¬ 
phragm motion, which are subjected to the boundary conditions 
z(a,0,t) = 0. Eliminating the time dependence e*"* from both sides 
of Eq. (17) leads to 


1 dz(r, 8) d 2 z(r, 9) 
7 dr + dr 2 


1 d 2 z(r, 8) 
+ r 2 dd 2 


+ K 2 z(r, 8) 


Ap P 

Pm C m 


(18) 


where m m , k m , £ m , F m and p m are the mass, stiffness, damping, ten¬ 
sion and surface density of the piezoelectric diaphragm. 

(a) m=0, n= 1 (b) m=0, n=2 




(c) m=l, n=l 


(d) m=l, n=2 






Fig. 4. The green’s function for a circular piezoelectric diaphragm fixed at the rim 
r = a, when a = 1.0, a = 1.0 and p = 0.3. 


Fig. 5. Measured dynamic response of the piezoelectric diaphragm before and after 
tuning: (a) resonant frequency shift and (b) spectrum of the generated voltage. 
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Table 1 

Piezoelectric transducer parameters and dimensions of the Rijke-Zhao tube. 


Weight (grams) 

Stiffness (N/m) 

Capacitance (nF) 

Resonant frequency (Hz) 

Diameter (mm) 

9.8 

5 x 10 3 

107 

290 

63.5 

a (degree) 

Length z u /z d 

Length z d3 /z d 

Length z d4 /z d 

Diameter A (m 2 ) 

15 

0.5 

2.41 

4.34 

0.0079 




frT' 


o. 

£ 






a. 

£ 




£ 





a, 

£ 





a. 

£ 



H 1 ' 


O. 

£ 



Time (s) 


Frequency (Hz) 


Fig. 7. Measured thermoacoustic oscillations and temperatures along the bifurcating branches of the Rijke-Zhao tube. 
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where c 2 m = r m /p m and k 2 = ^ The diaphragm damp¬ 

ing is generally sufficiently weak and so 1 m « 0. 

To solve Eq. (18), Green’s function approach is applicable and 
we can solve the partial differential equation 


\dg ctg \c?g 2 5(r -11)5(8 - 4>) 

r dr dr 2 r 2 dd 2 ^ r 


(19) 


where 0 < r,lZ < a and 0 ^ 0, </> ^ 2n, with the boundary conditions 
lim|g(r, 0\TZ, (j>)\ < oo g(a, 0\U, <t>) = 0 0^6,(/>^2n (20) 

r-> 0 


Since <5(0 - tj>) can be expressed in Fourier cosine series [41] as 


-i oo 

S{6 -<!>) = 2^J2 6m cos [ m ( e ~ <M1 

L m= 0 

It is anticipated that 

oo 

g(r, 6\n, <f>) = J2 e ™gm( r \K) cos (m(8 - </>)) 

m= 0 

where e m _ 0 = 1 and e m ^ 0 = 2. 

Substituting Eq. (22) into Eq. (19) yields 

1 9g m (7g m ( 2 m 2 \ S(r - TZ) 

r dr dr 2 \ r 2 ) Sm 2nr 


( 21 ) 


( 22 ) 


(23) 


The right hand side of Eq. (23) can be expanded as a Fourier-Bessel 
series as 


S(r - U) 
2nr 


1 

7ta 2 


y^ J m (k mr ,K/a)J m (k mn r/a) 

n=l Jm(kmn) 


(24) 


where k mn is the nth root of Bessel functions of order m of the first 
kind J m (k ) = 0. The prime denotes spatial derivative. Because the 
delta function is a Fourier-Bessel series, the solution to Eq. (23) is 
anticipated in the form of 

°° //c 

= (25) 

n=1 ^ ' 


Substituting Eqs. (25) and (24) into Eq. (23) gives rise to 

A . . 1 

*^mn — _ _7 / n \ 

710 (/<L/a 2 - ^ 2 )7m(fc m n) 

Therefore, the green’s function of the diaphragm motion is given as 


g(r,o\n, 4>) 


1 Jm(kmn'R'/ a ) 

(l7 mn /a 2 - K 2 )j*(k mn ) 


cos[m(0 - (j >)] 


(27) 


Fig. 4 illustrates the Green function g (normalized with its max¬ 
imum value) for a circular piezoelectric diaphragm fixed at the rim 
(i.e. r = a), as 0 - (j> is varied from 0 to 2n. It can be seen that there 
are multiple resonant modes of the piezoelectric diaphragm. And 
(a) m = 0,n = 1 corresponds the fundamental mode with lowest fre¬ 
quency («290 Hz). If the frequency of the oncoming pressure 
waves is close to one of the piezoelectric diaphragm resonant 
modes, then the diaphragm motion is maximized and so the volt- 
age/power output. 

Fig. 5 shows the measured dynamic response of the piezoelec¬ 
tric diaphragm before and after tuning. It can be seen from 
Fig. 5b that the piezoelectric generator can increase the voltage 
output by 105% in comparison with that without tuning. This is 
due to the fact that the tuning results in the shift of the resonance 
frequency of the diaphragm from 290 Hz to around 180 Hz, as 
shown in Fig. 5a. The spectrum indicates that in order to maximize 
the power output from piezoelectric generator, the thermoacoustic 
oscillations must be at frequencies closed to the resonant frequen¬ 
cies of the piezo-diaphragm vibration. 



Fig. 8. Temperature measurement by using an infrared camera. 




Fig. 9. (a) measured pressure mode shape along the upper longer branch and the 
mother tube and (b) Measurement and prediction of velocity oscillation near the 
piezo-diaphragm by using a hot wire anemometry. 


3. Description of experiment 

A Rijke-Zhao tube with a unique structure and geometry is 
designed and tested to demonstrate the thermoacoustic energy 
conversion. It has a mother tube (bottom stem), which splits into 
two daughter tubes (i.e. upper branches) with different lengths, 
see Fig. 6. As a heat source (flame, solar or electrical heater, etc.) 
is placed inside the mother tube, it provides a mechanism to pro¬ 
duce convection-driven thermoacoustic oscillation. The geometric 
and dimensions of the Rijke-Zhao tube and piezo diaphragm is 
shown in the Table 1. 
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Fig. 10. Measured power output of the piezo-diaphragm and the acoustical energy conversion efficiency with varying resistance R L in Ohms. 


The experimental setup is shown in Fig. 6. T-shaped junctions 1 
protruding from both bifurcating branches are designed for placing 
microphones and thermocouples. To measure pressure fluctuations 
along the Rijke-Zhao tube, two arrays of B&K microphones (Type 
4794) are implemented: one array with 4 microphones applied to 
the longer upper branch, the other with 3 microphones attached to 
the shorter branch. The measurements are shown in Fig. 7. It can 
be seen that there are thermoacoustic oscillations in both branches, 
and the maximum pressure oscillation can reach about 139 dB. The 
dominant mode in Rijke-Zhao tube is at around 176 Hz and it has 
several harmonics as shown in the bottom graph. This frequency cor¬ 
responds to the fundamental mode with a wavelength of twice of the 
total length of the bottom stem and the upper longer branch. 

The temperature of the acoustic fields in the bifurcating upper 
branches are measured by using two arrays of K-type thermocou¬ 
ples placed side by side to the microphones. This is denoted by 
the dash line in Fig. 7. It can be seen that one branch is associated 
with self-sustained ‘hot’ oscillations (approximately 450 K), while 
the other one is with ‘cold’ oscillations at ambient temperature (« 
300 K). This is due to the fact that ‘hot’ air goes through the upper 
longer branch, while the ambient air is ‘sucked’ into the other 
branch [39]. To further validate the temperature measurement, an 
infrared camera is used. It captures the quartz tube surface temper¬ 
ature spectrum. However, it can represent the ‘hot’ status of the 
sound waves inside each branch. The camera is set to operating 
range of 0-500 °C (corresponding to 273-773 K). Fig. 8 shows the 
measured temperature contour. It can be seen that the surface tem¬ 
perature of the short branch is around 26 °C («299 K), indicating 
that sound waves inside the branch are very close to room temper¬ 
ature. And this measurement is consistent with our thermocouples 
measurements. Note that the measurement error is ±2°. 


1 Note that the presence of T-shaped tubes might attenuate the thermoacoustic 
oscillations present due to boundary losses. However, they are not affecting the 

power output measurement from the piezoelectric membrane, since it is actuated by 
the pressure fluctuation at the end of the tube in the presence of T-shaped tubes. In 

practice, power output from the piezoelectric transducer might be increased by 

removing the T-shaped tubes. 


Finally, the ‘oscillating’ nature of thermoacoustics in the ‘cold’ 
branch is confirmed via direct measurement by using hot wire ane- 
mometry (HWA), as shown in Fig. 9b. The waveform deformation is 
most likely due to the presence of 2nd harmonic mode at about 
360 Hz, as revealed by conducting spectrum analysis to the veloc¬ 
ity measurement. The maximum volume flow rate is approxi¬ 
mately 4.77 x 10~ 4 m 3 /s. To gain insights on the acoustic field, 
mode shape along the longer upper branch is estimated and shown 
in Fig. 9a. The total length of the upper branch and the mother tube 
equals half of the wavelength of the thermoacoustic oscillation. In 
addition, good agreement between numerical prediction and the 
measured one is observed. This confirms that Rijke-Zhao tube is 
a reliable device to produce self-sustained thermoacoustic oscilla¬ 
tions to be utilized for energy harvesting. 


4. Experimental measurement of output power and efficiency 

The thermoacoustic oscillations at ambient temperature in the 
bifurcating branch enable the application of a piezoelectric gener¬ 
ator to harvest thermal energy. In our experiment, a piezoelectric 
diaphragm (Piezo Systems Inc. T216-A4NO-573X) made of Lead- 
Zirconate-Titanate is used. It has diameter of 63.5 mm, thickness 
0.41 mm and capacitance of 107 nF at resonant frequency of 
290 Hz. Before attaching the piezoelectric diaphragm to the end 
of the upper short branch as show in Fig. 6, it is tuned to maximize 
its electricity output by adding an aluminum disk-plate weighing 
2.85 g at the center of the diaphragm. 

The performance of the convection-driven thermo-acoustic- 
piezo system is shown in Fig. 10. It is can be seen from Fig. 10a that 
the maximum output power is about 2.1 mW. Compared with 
1.28 mW from the conduction-driven thermo-acoustic-piezo sys¬ 
tem proposed by Smoker et al. [23], our present system produced 
60% more power. In addition, tuning the piezo-diaphragm can 
greatly enhance the energy harvesting performance. Note that 
the electrical power output E L is calculated by using 

El = V 2 m JR L = V 2 /2 R l (28) 
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Fig. 11. A red LED is powered by the Rijke-Zhao thermo-acoustic-piezo system, which involves no heat exchangers and stacks (or regenerators). (For interpretation of the 
references to colour in this figure legend, the reader is referred to the web version of this article.) 


where V rms = V/y/2 denotes the root mean square of the output 
voltage V. The acoustical energy conversion efficiency, ;/ c , is shown 
in Fig. 10b. The efficiency is defined as r\ c = E L /p ■ V, where p and V 
denote the pressure and volume velocity at the piezo-diaphragm 
position. It can be seen that the tuned piezeoelectric diaphragm 
gives rise to the maximum energy conversion efficiency of approx¬ 
imately 22% when the load resistance is 5070 Q. Compared with the 
maximum efficiency 9.8% obtained by Smoker et al. [23], our ther¬ 
moacoustic system works 105% more efficient than the conduction- 
driven thermo-acoustic-piezo system. In addition, the overall 
efficiency of converting input heat power into output electrical 
power of conduction-driven system as designed by Smoker et al. 
[23] was shown to be about 0.00028%. In comparison with our pres¬ 
ent system, it is approximately 65% lower. The improvement of the 
overall efficiency is most likely due to the larger pressure oscilla¬ 
tions used to excite the piezoelectric generator and increased 
acoustical energy conversion efficiency. 

A practical demonstration is then made by using the electricity 
generated from the piezoelectric diaphragm to power a LED, as 
shown in Fig. 11. The successful application indicates that the 
Rijke-Zhao tube cannot only generate dual temperature thermoa¬ 
coustic oscillations but also provide a platform for energy harvest¬ 
ing. It is worth noting that our present energy conversion system 
involves no heat exchangers and stack (or regenerators). 


5. Discussion and conclusions 

In summary, a nonlinear thermo-acoustic-piezo model is firstly 
developed to gain insight on the energy conversion process in our 
Rijke-Zhao system. It captures two physical processes: the gener¬ 
ation and propagation of sound waves by inputting thermal energy 


and the response of the piezoelectric diaphragm to sound waves. 
Then, a convection-driven thermo-acoustic-piezo system for en¬ 
ergy harvesting is designed and experimentally tested. Unlike con¬ 
ventional conduction-driven thermoacoustic system, our present 
system is a convection-driven one. Furthermore, it does not involve 
any heat exchanger and stack (or regenerators). Comparing our 
present system performance with a similiar but conduction-driven 
thermo-acoustic-piezo system reveals that the power output is in¬ 
creased by 60%. And acoustical energy conversion efficiency is in¬ 
creased by 105%. The experimental results confirm that our 
system is much simpler in design, lower cost in fabrication and en¬ 
ergy-efficient. 
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Appendix A. Matrices U and V 

The coefficient matrices U and V are given as 

/ 0 0 1-1 0 \ 

0 0 0 1 -1 

0 0 At-ntt) -*| A.1I+M4) 

u= C2 c3 c4 

1+S?1 tt -l+M,(2-|-c D )-M[(l-§l-f) 0 0 0 

0 0 0 t 
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Appendix B. Mechanical-electrical coupling piezoelectric 
generator 


To maximize electric energy output, the piezoelectric generator 
operates at the resonance frequency of the vibrational structure to 
produce the maximum deformation of the piezoelectric materials. 
Near the resonance frequency, the mechanical behavior of the pie¬ 
zoelectric structure is well described by a one degree of freedom 
(1DOF) system. For simplicity, the single modal vibration of the 
mechanical structure is used to investigate the performance of 
the system. For a 1DOF system generator, such as a piezoelectric 
diaphragm vibrating under single bending modal, the mechani¬ 
cal-electrical coupling equations can be expressed as follows 


d z dz , 

m d? + C dt " kmeV = ApA p (29a) 

q - k me z ~^v=0 (29b) 

c p 


where m, c, k are modal mechanical mass, damping coefficient and 
stiffness respectively. k me is the modal piezoelectric coupling stiff¬ 
ness, A p is the modal external mechanical pressure applied, z is 
the modal displacement, v is the difference in electric potential be¬ 
tween the electrodes, q is the electric charge on the electrodes, and 
c p is the material capacitance. Taking Laplace transform to Eqs. 
(29a) and (29b) and manipulating them leads to the transfer func¬ 
tion of v(s) /Ap(s) as 


v{s) 

Ap(s) 


1 C km 


(30) 


where C(s ) = q(s)/v(s). It can be seen that due to the coupling be¬ 
tween the mechanical and electrical components, the exciting force 
give rise to the displacement in the piezo, which induces the elec¬ 
trical charge on the electrodes. This explains how the electrical 
power is generated from piezoelectric diaphragm as we used in 
our experiment. 
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